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Abstract. We investigate for the first time the potential of angular auto- and cross-correlation 
power spectra in identifying sterile neutrino dark matter in the cosmic X-ray background. 
We take as reference the performance of the soon-to-be-launched eROSITA satellite. The 
main astrophysical background sources against sterile neutrino decays are active galactic 
nuclei, galaxies powered by X-ray binaries, and clusters of galaxies. While sterile neutrino 
decays are always subdominant in the auto-correlation power spectra, they can be efficiently 
enhanced when cross-correlating with tracers of the dark matter distribution such as galaxies 
in the 2MASS catalogues. We show that the planned four-years eROSITA all-sky survey will 
provide a large enough photon statistics to potentially yield very stringent constraints on the 
decay lifetime, enabling to firmly test the recently claimed 3.56-keV X-ray line found towards 
several clusters and galaxies and its decaying dark matter interpretation. However, we also 
show that in order to fully exploit the potential of eROSITA for dark matter searches, it is 
vital to overcome the shot-noise limitations inherent to galaxy catalogues as tracers for the 
dark matter distribution. 
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1 Introduction 


Revealing the non-gravitational nature of dark matter particles is one of the most impor¬ 
tant goals of modern astrophysics and cosmology. Even though weakly interacting massive 
particles (WIMPs) are among the most popular and well-motivated candidates [1, 2], so 
far no established signature thereof has been found at colliders, direct, or indirect searches. 
This urges the investigation of alternative possibilities. Sterile neutrinos with keV masses 
are well-known and well-motivated candidates for dark matter [3-6]. Observationally, keV- 
rnass sterile neutrinos are interesting as they behave as warm dark matter (WDM). This 
can potentially mitigate some of the small-scales problems of the cold dark matter (CDM) 
scenario, e.g., the so called missing satellite and too-big-to-fail problems, or alternatively 
provide indirect evidence for sterile neutrino dark matter (see, e.g., [7]). 

Sterile neutrinos are spin- 1/2 singlets under the Standard Model gauge group, and can 
decay into photon-neutrino pairs, u s —> 71 /. A clear smoking gun signal for sterile neutrino 
dark matter would be the detection of monochromatic photons at half of the dark matter 
mass. The corresponding decay rate is proportional to the mixing angle 0 2 between active 
and sterile neutrinos, and is given by [5, 8 ] 


T Va ~ (7.2 x 10 29 s )" 1 


/ sin 2 29 \ t m Us \ 5 
V 10 “ 8 ) VI keV/ ’ 


( 1 . 1 ) 


where m„ a is the sterile neutrino mass. 
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Monochromatic lines from sterile neutrinos have been searched for in various different 
targets, such as nearby galaxies and galaxy clusters [9-14]. Recently, there have been claims 
of the detection of an unidentified 3.56-keV line from a stacked sample of galaxy clusters [15], 
and from Andromeda and the Perseus cluster [16], with a subsequent number of works on the 
issue, some of which confirmed the claim [17-20] and some not [21-27]. This is an ongoing 
debate that deserves further investigation (see, e.g., [28]). 

One modern and very promising approach to indirect dark matter searches with photons 
at all wavelengths is the cross-correlations between the electromagnetic signal and tracers of 
the dark matter distribution [29, 30]. Up to now, these searches were mostly related to dark 
matter annihilation in the gamma-ray regime, with few studies also on decaying dark matter 
[29, 31, 32]. The main reason being that the focus of indirect dark matter searches so far has 
been mainly on WIMPs and their gamma-ray annihilation and decay products. However, as 
already mentioned, the lack of any strong evidence for WIMPs have recently broadened the 
interest of the community, particularly after the 3.56-keV line claim and its sterile neutrino 
dark matter interpretation. 

Cross-correlation searches aim at the dark matter signal coming from cosmological dis¬ 
tances, and have the advantage of exploiting simultaneously spatial and spectral information 
from the full sky. As tracers for dark matter in the local Universe, galaxy catalogues [30] and 
cosmic shear (weak gravitational lensing) [29, 33] proved to be very promising. Addition¬ 
ally, with a tomographic approach one can further exploit the redshift information, which 
allows a better discrimination between various astrophysical and dark matter related photon 
sources [31, 34, 35]. The power of such cross-correlation searches manifested recently in the 
first discovery of cross-correlations of the Fermi -Large Area Telescope gamma-ray sky with 
the CMB [36] and galaxy catalogues [37], with important implications for WIMP dark matter 
searches studied in Refs. [38, 39]. A recent exhaustive review about the gamma-ray emission 
from cosmological distances, and the potential of cross-correlation studies, can be found in 
Ref. [40] (see also Ref. [41] for a detailed discussion about the expected WIMP dark matter 
signal). 

In the present paper, we investigate for the first time the potential of an angular power 
spectrum analysis of the cosmic X-ray background (CXB) in identifying sterile neutrino dark 
matter. The methods that we present are very general and can be used for all scenarios 
with decaying dark matter. We will show that the auto-correlation power spectrum of X-ray 
emission from sterile neutrinos is completely dominated by background sources, in particular 
by active galactic nuclei (AGNs) and clusters of galaxies. However, as we will demonstrate, 
the cross-correlation with tracers of the dark matter distribution in the Universe, like large 
galaxy catalogues, can greatly enhance the sensitivity towards a sterile neutrino contribution 
to the X-ray flux, and hence can be a sensitive probe of dark matter signals. 

For the purpose of such a CXB analysis, a deep and (nearly) all-sky survey with good 
energy and angular resolution is a fundamental requirement. The last X-ray all-sky survey 
has been performed by ROSAT (e.g., [42]). Current X-ray instruments such as XMM-Newton 
(e.g., [43]) and Chandra (e.g., [44]) have impressive performances and achievements. However, 
they only operate in pointed observations. The next X-ray all-sky survey will be performed by 
eROSITA [45], whose launch is scheduled for 2016, with a set of 6-month full-sky observations 
up to four years. Therefore, eROSITA will be the reference for our predictions. 

As an exemplary (and readily available) tracer for the nearby dark matter distribution 
in the local Universe, we will adopt the 2MASS galaxy catalogue [46]. It provides a nearly 
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complete information on the distribution of galaxies up to z ~ 0 . 1 , with a sky coverage of 
91%. We will discuss also the potential ultimate reach of our method when other tracers are 
adopted. 

The paper is organised as follows. In section 2 we will discuss in detail the different 
background and the signal contributions to the CXB. Their auto-correlation power spectrum 
is discussed in section 3, whereas the cross-correlation power spectrum with tracers of dark 
matter is discussed in section 4. In section 5 we present our projected sensitivity for future 
instruments. Eventually, in section 6 , we present our conclusions. 

If not specified otherwise, we adopt a halo mass Ma defined with respect to a density 
that is A = 200 times the critical density of the Universe at redshift z. Throughout the 
paper, we assume a standard ACDM cosmology with the parameters = 0.73, H m = 0.27, 
= 0.05, Qdm = 0.22, and Hq = 100 /ikm s _1 Mpc -1 with h = 0.7. 


2 Signal and backgrounds in X-ray searches for sterile neutrinos 

In the following three sections, we will discuss and characterise the mean intensity of a 
sterile neutrino dark matter signal and of the most relevant background components. These 
backgrounds are fluxes from AGNs, X-ray binaries hosted in galaxies, and emission from 
clusters of galaxies. 


2.1 Sterile neutrino dark matter 

The intensity of photons from sterile neutrino decay into 'yi/ final states along the direction 
n, defined as a number of photons received per unit area, time, solid angle, and energy range, 
is given by 1 


I Us (E, xn) 


dz 

J 0 H(z ) (1 + z) 3 


Pv 3 ( Z 1 X n ) < 5 d 


(1 + z)E 


rjW 

2 J ’ 


( 2 . 1 ) 


where 5d is the Dirac delta function, H(z) = Hoy/£l m (l + z) 3 + Da is the Hubble parameter 
as function of redshift z, Pu 3 (z, x n ) = Pd m (u x n ) denotes the mass density of sterile neutrino 
dark matter in the direction n at a distance z, and x is the co-moving distance. 

We will first evaluate the mean intensity of extragalactic sterile neutrino dark matter 
and hence the overall contribution to the CXB. We note that adopting the ensemble average 
yields (p Us (z,x n )) = (1 + z) 3 £l Vs p c = (1 + z) 3 f 2 Hm £>r, where p c is the present critical density 
of the Universe. Therefore, after taking the ensemble average of Eq. (2.1), and convolving 
with a normal distribution in energy in order to represent the finite energy resolution of an 
X-ray telescope, the mean intensity due to the sterile neutrino decay is given by 

I Us (E)= d X W„ s ([l + z\E,z) , (2.2) 

Jo 

where the window function W Us (E, z) is 


W v .{E,z) 


_ ^dmPcr l/ s 

2(2TT) 3 / 2 m Us (l + z)(Je P 


(E-m v J 2 ) 2 ' 

2(1 + zYg I . 


(2.3) 


1 The lower limit for the redshift integration in Eqs. (2.1) is chosen to be 0 = 0.003. The mean intensity, as 

well as the angular power spectrum at high multipoles, are shown to be insensitive to this choice (e.g., [47]). 

See also the window function dependence on z in the right panel of Fig. 3. 
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Figure 1. Contributions to the CXB mean intensity as a function of energy. We show the sterile 
neutrino signal for our reference scenario as described in the text (motivated by the findings of 
Refs. [15, 16]), compared with other astrophysical contributions: unresolved AGNs and galaxies, and 
galaxy clusters (both resolved and unresolved). 

We adopt ge = 0.138 keV/\/81n2 = 0.0586 keV corresponding to the energy resolution 
expected for eROSITA [45] at 6 keV. 2 As expected, Eqs. (2.2) and (2.3) show that, for a given 
energy E, the dominant contribution come mainly from redshifts around z = m Vs /2E — 1. 

Fig. 1 shows the mean intensity due to sterile neutrino decays, for a mass 7.12 keV 
and sin 2 26 = 7.6 x 10~ 11 , which will be our reference scenario motivated by the findings 
in Refs. [15, 16], compared with the CXB contribution from AGNs, galaxies, and galaxy 
clusters (these background components will be discussed in sections 2.2-2.4). The sterile 
neutrino component is completely subdominant, being smaller by two orders of magnitude 
with respect to the dominant AGN contribution, even at the peak energy around 3.5 keV. 

2.2 Active galactic nuclei 

AGNs are believed to provide the dominant contribution to the measured CXB [48-50]. 
We are here interested in the contribute from unresolved AGNs, since we can assume that 
all resolved point sources can be masked before performing the angular power-spectrum 
analysis. The unresolved contribution is extrapolated from the resolved one. Deep X-ray 
surveys suggest that the unresolved component contributes 76% to 82% in the 0.5-2 keV and 
2-8 keV energy bands to the overall CXB, respectively [51]. 

We describe AGNs via their X-ray luminosity function (XLF) defined as the co-moving 
number density per unit logarithmic luminosity range 4 >agn(Fx,- 2) = dttAGN/dlog 10 Lx- 
There are several parameterisations for the XLF and we adopt the Luminosity And Density 

2 Value taken from http://www.mpe.mpg.de/455799/instrument. We approximate the energy resolution 
as being constant in our energy range of interest, while in reality it improves at lower energies and worsen 
at higher energies. For example, in the case of XMM-Newton (see the XMM-Newton Users’ Handbook), the 
FWHM energy resolution at 1, 6 and 10 keV is about 70, 150 and 180 eV, respectively. 
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Evolution (LADE) model [52]: 


72 "I -i 


®agn(Lx,z) = K{z) 





(2.4) 


where K(z) = io _4 - 53_0 - 19 ( 1+z ) Mpc^ 3 , 71 = 0.62, and 72 = 3.01. The characteristic lumi¬ 
nosity L it (z), where the power-law break happens, depends on redshift and is parameterised 
as 


L*{z) = Lq 


f 1 + z c 

v 1 + z 


/1 + z c 



(2.5) 


where Lq = 10 44 7 ' erg s _1 , z c = 0.75, p\ = 6.36, and P 2 = —0.24. Note that this XLF refers 
to luminosities integrated between E m i n = 2 keV and E max = 10 keV in the AGN rest frame. 

The purpose of the current study is not to get a precise estimate of the mean intensity 
of unresolved AGNs, but to show the power of angular fluctuation to identify sterile neutrino 
dark matter. Therefore, we simply assume that the each AGN features a simple power-law 
spectrum, E ~ r , with a spectral index T = 1.7 in the relevant X-ray energy range (e.g., [53]). 
With this assumption, the differential luminosity at an energy E (defined as a number of 
photon emitted per unit time and per unit energy range), at redshift z, is 


£x(E, z) = L x 2 _ r 2 _ r E r , (2.6) 

^max ^max 

where Lx denotes the luminosity used in Eq. (2.4). The mean AGN intensity is then calcu¬ 
lated, similarly to Eq. (2.2), as 3 


Iagn{E) = / dx Wagn([1 + z\E, z) , 

Jo 

where the AGN window function is given as 

1 f LX ’max 


Wagk(E, z) = 


4vr In 10 7 Lx , 


Lx 


^AGN (Lx, z)Cx(E, z) 


(2.7) 


( 2 . 8 ) 


For the lower and upper limit of luminosity integration, we adopt Tx,min = 10 4i erg s 1 and 


Lx, r 


Z = 


4vrci j E S ens 


E 2 ~ r 

max 


- E, 


2 -r 


(1 + z?-r E^r - El^ 


(2.9) 


respectively (e.g., [52]). F sens is the flux sensitivity of the considered X-ray telescope, and we 
adopt 4.4 x 10 ~ 14 erg cm -2 s ^ 1 for the energy range between Emin = 0.5 keV and E max = 
2 keV for eROSITA [45], which corresponds to the point-like sensitivity expected after the 
first six months of survey . 4 The mean intensity of unresolved AGNs computed in this way is 
shown in Fig. 1. It is by far the dominant contribution to the CXB. 

The AGN luminosity density of the adopted model has a 1 cr uncertainty below 50% at 
low, z < 0.2, and high, z > 3 redshifts, and below 20% in the middle range [52], There¬ 
fore, while, as explained above, we refrain from a precise estimation of the unresolved AGN 
contribution, we note that our conclusions are robust against the present uncertainty. 

3 The lower limit for the redshift integration in Eq. (2.7) for AGNs, as well as for galaxies in section 2.3, 
is chosen to be 10~ 4 . This does not impact the calculations at all because all the AGNs and galaxies within 
this redshift give fluxes above the eROSITA sensitivity and hence they are completely resolved. See also the 
window function dependence on z in the right panel of Fig. 3. 

4 This implies that we will obtain slightly more conservative results than considering the sensitivity for the 
full four year survey, but we tested that the change in our final constraints is not significant (mostly since the 
result is limited by other factors, see discussion below). 
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2.3 Galaxies 


Galaxies are also X-ray sources as they host X-ray binaries. The flux distribution of galaxies 
approaches that of AGNs at the faint end [51]. We adopt the XLF (in the rest-frame energy 
range band 0.5-2 keV) from Ref. [54], where the redshift evolution was studied from a galaxy 
sample below z ~ 1. The galaxy XLF is well represented with the following log-normal 
distribution function: 


$gai {Lx,z) = < 3 ?* 



1—a 

exp 




( 2 . 10 ) 


where = 10~ 2 ' 23 Mpc -3 , a = 1.43, and a = 0.72. The characteristic luminosity L* 
evolves with redshifts while 4>* is constant (pure-luminosity evolution), and its evolution is 
parameterised as L* = 10 39 ' 74 erg s~ 4 [(l + z)/1.25] 1-9 . 

The CXB mean intensity I sa \(E), as well as its window function W ga \(E, z), are com¬ 
puted in the same way as for the case of AGNs, by using the galaxy XLF 4> ga i(Lx) z) in 
Eqs. (2.7) and (2.8). This time, however, we adopt a E~ 2 spectrum [55], and thus the 
differential luminosity is given by 


£x(E,z) = Lx 


E 


-2 


ln(F/ max /Emin) 


( 2 . 11 ) 


where E m i n = 0.5 keV and E max = 2 keV, and the upper limit of the luminosity integration 
is simply computed as Lx ma (z) = 47rd|F sen s, with F S P = 4.4 x 1CU 14 erg cm -2 s~ 4 as for 
AGNs. The mean intensity of unresolved galaxies computed in this way is shown in Fig. 1. 
It is only about 10% of the emission from unresolved AGNs, but still an order of magnitude 
above our reference sterile neutrino signal. Also in this case, the 1 a uncertainty on the 
luminosity density is below 50% [54] and, therefore, our conclusions are robust. 


2.4 Clusters of galaxies 

Galaxy clusters are prominent X-ray emitters as the ambient gas, the intra-cluster medium 
(ICM), is a hot thermal plasma radiating bremsstrahlung emission (see [56] for a review). 
Indeed, X-ray observations are one of the main method to identify clusters. Clusters also 
represent some of the most promising targets to identify sterile neutrino dark matter, or any 
decaying dark matter candidate, as this process is directly proportional to the halo mass and 
they indeed are the most massive halos in the Universe. 

Future X-ray instruments will help in narrowing down the contribute of galaxy clusters 
to the CXB. In fact, eROSITA is predicted to be able to resolve all galaxy clusters in the 
Universe [45]. However, differently from the previous cases of AGNs and galaxies, we must 
here consider the whole population of resolved and unresolved clusters and their thermal 
X-ray emission. This is because one of the questions here is: do clusters of galaxies also 
represent the bulk of the angular power spectrum for the sterile neutrino dark matter, and, 
therefore, is their thermal bremsstrahlung emission to be considered a dominant background 
in this approach? We underline that the bremsstrahlung emission depends on the gas density 
squared, which implies that significant morphological differences between a dark matter decay 
signal and the gas emission are expected. 


- 6 - 







Following the above notation, the mean cluster X-ray intensity due to bremsstrahlung 
radiation is given by 5 

r oo 

Id(E) = / dx W d ([l +z\E,z)(p 2 gas ) , (2.12) 

Jo 

where p gas is the ICM density, and 

<p «“ > = (s^) I dM200 Jito f dVp t “ {r|M ™)' (213) 


with integration boundaries that are discussed below, 
function is 


W cl (E, z) 


(fjfe Pc ) 2 

47T 


k s 


(k B T gas )- 1/2 

E 


The cluster bremsstrahlung window 

exp ( _ %^;)' (214) 


where k B T geLS is the ICM temperature in keV. We define kg = 9.5 x 10 32 gg- (adapted from [58]), 
where gs = 1.1 is the Gaunt factor for the free-free emission, such that I c i(E) is expressed 
in cm~ 2 s -1 sr -1 keV -1 . 

Our expression for the bremsstrahlung emission do not consider atomic line transitions 
in the ICM. This simplification is needed in order to treat the problem analytically. Atomic 
line transitions depend on temperature and mass of the cluster, and on the density squared 
as bremsstrahlung. The exclusion of atomic lines will hence have no impact on our study 
of morphological differences in the angular power spectrum for the different components. 
One might worry that our side-band analysis of section 5, which takes into account spectral 
information, is significantly affected by this approximation. However, as we will see, even 
neglecting entirely galaxy clusters in the analysis does not change much the final results. 
Hence, we believe that avoiding the transition lines does not impact much on this first study 
that we attempt here. 

The adopted halo mass function dn/dM 200 in Eq. (2.13), which will be later also used 
for the sterile neutrino angular power spectrum calculation, is based on the Tinker formalism 
[59] and obtained through the online application from [60]. In Eq. (2.13), the mass function 
is integrated above x h = 10 14 Mq, while the volume integral is within i? 200 - The 

transition between galaxy groups and galaxy clusters is not sharp or well-defined, and with 
the chosen we already include large galaxy groups in the integration. 

Finally, we need a prescription for the ICM density and temperature. We adopt the 
phenomenological model of [61], which is based on X-rays observations, and allows to assign 
a gas density (and temperature) to any galaxy cluster using its mass only, in such a way that 
the observed X-ray and Sunyaev-Zel’dovich scaling relations are correctly reproduced within 
about 20%. Following [61], we also statistically divide the cluster population in half merging 
clusters and half relaxed clusters, which are characterised by centrally cored and centrally 
peaked gas densities, respectively. This prescription needs however the clusters’ mass as M500 
which we obtain from M 200 using the method in [62]. As a result, we find that the mean 
intensity of galaxy clusters, as shown in Fig. 1, is at the level of unresolved galaxies. 

s The lower limit for the redshift integration in Eq. (2.12) is chosen to be z = 0.01 as the closest known 
galaxy clusters (see, e.g., [57]). 
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3 Auto-correlations of the angular power spectrum 


Discussing the auto-correlation angular power spectrum of the different signal and back¬ 
ground contributions to the CXB is useful for identifying the components that will yield the 
strongest contrast at a given angular scale 8 = it/ l, where t is the multipole. As usual, we 
define the intensity angular power spectrum of the signal and background fluxes as 

Q = (|«£m| 2 ) J (3-1) 


where the a£ m are given by the decomposition of the flux into spherical harmonics Ye m , 
namely 

atm = J dQ n I{n) Y/ m (n) . (3.2) 

The auto-correlation angular power spectrum of a given source population A (= n s , AGN, 
gal, cl) is then computed as [34, 63] 

Cf(E) = r% W A ([ 1 + z]E, z) 2 P A (k = -,zY (3.3) 

J o X V X J 


where Pj^{k^z) is the power spectrum of the source A at wave number k and redshift z. In 
the halo model, the latter can be divided into one- and two-halo terms, P\ = P^ h + P'f^' 
(higher halo terms are negligible). In the case of dark matter decay, it is the nonlinear matter 
power spectrum (P Us = P$) modelled as [64] 


P, lh = 


P, 2h = 


1 


1 


dM 2 00 
dM 2 oo 


dn 


dM .200 
dn 


4irr 2 drp dm (r) 


sin(/cr) 

kr 


1 2 


(3.4) 


-jjj—b(M 200 ,z) j 4:irr 2 drp dm (r) Sm ^ 


Phn{k, x) j 


where the mass integration starts from M^qq 1 '™ x h = 10 6 M 0 , the radial integration goes up 
to i?200) P\m(k,x) is the linear matter power spectrum (obtained from [60]), and b(M 2 Q dl z) 
is the linear bias [65]. The lower mass limit for the integration of the sterile neutrinos, 
which represent effectively a WDM candidate, can be chosen to be around 10 6 -10 8 M 0 
(e.g., [66, 67]). Since the precise choice of the lower mass limit does not affect the auto- and 
cross-power spectra that we will discuss, we adopt for definiteness ^Q lim X h = 10 6 M 0 . 
Note also that our Pu n (k,x) is obtained for a CDM scenario, but as the WDM halo mass 
function starts to drop below its CDM counterpart around 10 9 M 0 [67], this choice of P]j n 
has no impact on our conclusions. The negligible relevance of the two above choices, M^ 110 
and Pn n (k,x), is made clear by Fig. 7 which shows that the dominant contribution to the 
sterile neutrino power spectrum comes from halos with M200 X h > 10 11 M 0 . 

In order to calculate the angular power spectrum of Eqs. (3.5), we need the distribution 
of dark matter in halos p dm . For the present study, we will adopt the Navarro-Frenk-White 
(NFW) profile [68], which is given by 


Pdm(r ) 


_ Ps _ 

(r/r s ) (r/r s + l) 2 ’ 


(3.5) 


where p s and r s denote the scale density and radius, respectively. We use throughout the 
parameters M200 and P200 to indicate the halo mass and radius, respectively. Hence, our 
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Figure 2. Bias of the X-ray AGNs and galaxies as a function of redshift 2. 


scale radius is defined as r s = R 200 /C 200 , where 0200(^200, z) is the concentration parameter 
[69]. With this, one can show that, for a given halo mass M 200 , the scale density is given by 


_ M200 
47rrf 


ln(l + c 2 oo) 


C200 

1 + C200 



(3.6) 


In the case of clusters of galaxies, Pc\(k,x ) is similar to that of Eqs. (3.5), but with 
AC m X h = 10 14 Mq. Additionally, /9dm( r ) has to be substituted by p gSlS (r) 2 and (l/fidmPc) 
by (1/fibPc) 2 , and we assume that Pii n (A:, x) is related to the cluster power spectrum via the 
linear bias. 

We assume that AGNs and galaxies are good tracers of the dark matter density, and 
thus -PAGN,gai(&, z) = ^agn gal ( z )Rs{k, z), where we again use Eq. (3.5) for the nonlinear 
power spectrum of dark matter density fluctuation 5. For the AGNs, we adopted the halo 
linear bias from [65] and assume that the X-ray AGNs reside in dark matter halos with mass 
of 10 131 Mq [70]. For the X-ray emitting galaxies, we adopt the prescription given in [71]. 
The bias models for both AGNs and galaxies are shown in Fig. 2 as a function of redshift. 
The uncertainty on the AGN bias can be estimated to be below about 10%, while for the 
halo mass in which X-ray AGNs are supposed to reside it is of about 30% [70]. The model 
for the bias of the X-ray galaxies of Ref. [71] well reproduces measurements of star-forming 
and high-z Ly-break galaxies which in turn have uncertainties below about 20%. Note that 
the galaxy bias is very close to 1 for 2 < 1, which is the redshift range most interesting for 
dark matter searches. 

Finally, since AGNs and galaxies are point-like sources, there is an additional shot-noise 
contribution to the angular power spectrum, which is independent of angular scale t. This 
Poisson term is computed by (e.g., [63]) 


CpGN.gal^) 


1 [°° dx 

(4tt ) 2 In 10 J 0 x 2 



dLx 

Lx 


$ AGN,gal (Ax, z)Cx{E, Z ) 2 . 


(3.7) 


So far, we discussed only the differential (in energy) angular power spectrum, but the 
same argument applies to the angular power spectrum integrated over a certain energy range. 











Figure 3. Left. Auto-correlation angular power spectrum of CXB due to sterile neutrino decays for 
our reference model, (unresolved) AGNs, (unresolved) galaxies, and (resolved and unresolved) clusters 
of galaxies in the 3.4 —3.6 keV energy band. For AGNs and galaxies, both the correlation term (solid) 
and total including Poisson term (dotted) are shown. Right. Redshift-weighted window functions 
integrated in the 3.4 — 3.6 keV energy band and normalised at their maximum for sterile neutrino 
decays with our reference model, AGNs, galaxies, clusters of galaxies (fixing the ICM temperature to 
5 keV), and 2MRS and 2MXSC catalogues. 


One first integrates the window function Wa over the relevant energies, and then computes 
the power spectrum using Eq. (3.3). Similarly, for the Poisson term, one integrates £x over 
the chosen energy range and then computes the shot noise with Eq. (3.7). In order to em¬ 
phasise the difference, we will in the following explicitly show the energy dependence, Cg(E), 
when refereeing to differential power spectra, while the energy dependence is suppressed, Cg, 
for the integrated quantities. 

In the left panel of Fig. 3, we show the auto-correlation angular power spectrum inte¬ 
grated over the 3.4-3.6 keV energy band, adopting our reference model for the sterile neutrino 
decay. The component coming from sterile neutrinos is completely subdominant with respect 
to contributions from clusters at all multipoles, as well as from AGNs, and, to a lesser extent, 
from galaxies, for most angular ranges. Furthermore, we also show that at larger multipoles 
the contributions from the Poisson term become increasingly important and start to domi¬ 
nate. In the right panel of Fig. 3, we show the shape of the window function integrated in 
the 3.4-3.6 keV energy band for the different components, and for the 2MRS and 2MXSC 
galaxy catalogues that will be introduced in the next section. 

Actually, the power spectrum of the sterile neutrino component can be obtained by using 
a non-linear matter power spectrum directly calculated with, e.g., the Code for Anisotropies 
in the Microwave Background (CAMB; http://camb.info/). We adopted the halo model ap¬ 
proach in order to have control on the spatial part of the angular power spectrum calculation 
for clusters of galaxies, and kept this approach in all cases for consistency. In Fig. 4, we com¬ 
pare the sterile neutrino auto-correlation angular power spectrum obtained with the halo 
model approach to that obtained with the non-linear matter power spectrum from CAMB 
and find reasonably good agreement, to within a factor of less than two. 
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Figure 4. The auto-correlation angular power spectrum of sterile neutrino decays computed with the 
halo model (used in our main analysis) and with the non-linear matter power spectrum from CAMB, 
respectively. 


4 Cross correlations with the 2MASS galaxy catalogues 


We will now show that, since the signal from sterile neutrino decays follows the distribution 
of dark matter in the Universe, a cross-correlation between the X-ray signal and tracers of the 
dark matter distribution is a very promising strategy to isolate the dark matter contribution. 
To this end, we will compute the cross-correlation with large galaxy catalogues from the 
2MASS survey [46], which in fact trace the dark matter distribution in the local Universe. 

The 2MASS catalogue provides nearly complete information on the galaxy distribution 
for the local Universe. The 2MASS Redshift Survey (2MRS) [46] is based on spectroscopic 
redshift determination of -Admits ~ 43500 galaxies up to z ~ 0.1 with a sky coverage of 91%. 
The 2MRS galaxies follow approximately the redshift distribution given by 


dN 2MR S 
dz 


oc z exp 


Z2MRS 


(4.1) 


where Z2MRS = 0.033. For comparison, we also adopt the 2MASS Extended Source Catalogue 
(2MXSC), an even larger sample that contains IV2MXSC ~ 770000 galaxies with a sky coverage 
of 67% [72, 73], but with less accurate redshift determination through photometry. The 
redshift distribution of 2MXSC is approximately given by [73] 


^-A^mxsc 

dz 


1 9 

oc z ' exp 


Z2MXSC 


1.75'' 


(4.2) 


where £2MXSC = 0.07. We fix the normalisation of both the 2MRS and 2MXSC catalogues 
such that the redshift integration gives the total number galaxies A^mrs and Wmxsc , re¬ 
spectively. 

The angular cross-power spectrum of a given source population A (= is s , AGN, gal, cl) 
with the 2MRS or 2MXSC galaxies, labelled B, is therefore computed as [63] 

Cf’ B (E) = J ^W A ([l + z)E,z)W B ( X )PA,B(k = (4.3) 
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Figure 5. Cross-correlation angular power spectrum of CXB due to sterile neutrino decays for the 
reference model, (unresolved) AGNs, (unresolved) galaxies, and (resolved and unresolved) clusters 
of galaxies in the 3.4-3.6 keV energy band. For sterile neutrinos, we also show the case when the 
integration upper mass limit is fixed to x h = 10 13 M 0 (see main text for details). We overlay 

to the cluster component the error bars AQ/VAf, with A£ = 1 /2 the bin size, estimated from the 
diagonal part of the covariance matrix for the total power of Eq. 5.1. Left: Cross correlation with 
2MRS. Right: Cross correlation with 2MXSC. 


where VFb(x) = (dz / dx)(dNs / N^dz) is the galaxy catalogue window function (normalised to 
unity), and PA,B(k,z) is the cross-power spectrum. For the cross-power spectrum, we adopt 
for the galaxy catalogues bias &2MRS/2MXSC = b ga \. In case of cross-correlation with clusters 
of galaxies, the cross-power spectrum, P c i b = -P c Vb + P c 2h B , reads 
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P\in{k,x)- (4-5) 


In the case of cross correlation with sterile neutrino dark matter, the cross-power spec¬ 
trum reads as in Eqs. (4.4) and (4.5), but one has to substitute /Jg as (r) with pdm(i’) and 
(l/Ob/0c) 2 with (1/fldm/Oc)- F° r tl ie cross-correlation with X-ray AGNs or galaxies, we have 
PAGN/gal,B(P) z) = ^AGN/gal(*)^2MRS/2MXSC(^) P&{k, z'). 

We show in Fig. 5 the cross-correlation angular power spectrum, ’ , cross-correlated 
with 2MRS (left panel) and 2MXSC (right panel) galaxies, integrated over 3.4-3.6 keV, and 
compare the sterile neutrino component with AGNs, galaxies and clusters. The shapes of 
the power spectra are similar for all components, as they trace the dark matter distribution 
quite well, while some deviations are observed for clusters where the emission depends on 
the gas density squared. 
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Figure 6. Left: Auto-correlation power spectrum of sterile neutrino decays, for the reference model, 
and clusters of galaxies in the 3.4-3.6 keV energy band. For clusters of galaxies, we show the case 
when all resolved and unresolved objects are included, and the cases where we exclude the objects 
exceeding the eROSITA expected sensitivities (for point-like/extended sources and eRASSl/eRASS8 
results). For sterile neutrinos, we also show the case when the integration upper mass limit is fixed 
to M^qo™ x h = 10 13 Mg. Right: The same but for cross-correlation with 2MRS. 

We find that the cross-correlation signal of (resolved and unresolved) galaxy clusters 
dominates over all components at all multipole scales. At the same time, the sterile neutrino 
component, subdonrinant in the auto-correlation, is now comparable to the AGN and galaxy 
components. In particular, the sterile neutrino component dominates over the AGNs and 
galaxies over all multipoles for the cross correlation with 2MRS. When cross-correlating with 
2MXSC, AGNs and galaxies dominate over the sterile neutrino component at all multipoles 
above l ~ 5, while being still comparable. This can be explained by the fact that unresolved 
AGNs and galaxies grow in number with redshift, therefore correlating more with the higher 
redshift sample of 2MXSC. At the same time, most of the power from sterile neutrinos comes 
from low redshifts, which causes the 2MXSC cross-correlation power spectrum to be lower 
with respect to the 2MRS case. This can be seen also looking at the redshift variations of 
the window functions of different components in the right panel of Fig. 3. 

We also try an additional approach to further reduce the contribution from clusters. This 
is simply to exclude completely clusters of galaxies from such an analysis. One would think 
that by doing so, we would exclude most of the sterile neutrino signal too, but we will see that 
the contrary is true. In Fig. 6, we show how the contribution from clusters of galaxies can be 
dramatically lowered when excluding clusters that will be resolved by eROSITA. We consider 
both point-source and extended expected sensitivities [45] to bracket the case of extended 
and smaller, far-away point-like clusters, from the first six months of survey (eRASSl) and 
the complete survey (eRASS8). The left panel of Fig. 6 shows the auto-correlation power 
spectra, while the right panel shows the cross-correlation with 2MRS. We compare these with 
the sterile neutrino component, and with the same sterile neutrino component but integrated 
only up to x h = 10 13 M 0 , effectively excluding all clusters of galaxies. The figures 

clearly show that excluding clusters of galaxies from such an analysis does not have a big 
impact on the magnitude of the sterile neutrino component, and, at the same time, that 
with the eROSITA survey it will be realistic to exclude them as all will be resolved. The 
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relative importance of galaxy clusters becomes even clearer when looking at the individual 
contributions from different halo sizes in Fig. 7 where we show that a large contribution 
to the sterile neutrino power spectra comes from galaxy-size halos. For completeness, we 
show the sterile neutrino component integrated only up to g™ x h = 10 13 Mq also in the 
cross-correlation spectra of the above Fig. 5. 

To conclude this section, in Fig. 8, we show the 2MRS (left) and 2MXSC (right) cross 
correlation (E) at l = 50 as a function of energy. The cluster component dominates at all 
energies. The sterile neutrino component dominates over those of AGNs and galaxies around 
its peak energy for 2MRS, while it is just below for 2MXSC. This is due not only to the 
fact that high-redshift AGNs and galaxies cross-correlate more with the higher redshifts of 
2MXSC, as discussed for Fig. 5, but also to the fact that, in this case, the sterile neutrino 
decay line is broader due to these higher redshifts too. Looking back at the 2MXSC cross 
correlation in the right panel of Fig. 5, we can argue that a energy range broader than 
3.4-3.6 keV, better including the sterile neutrino decay line peak, would increase its cross¬ 
correlation with 2MXSC, but would also increase that of the other components. Note that 
the choice of t = 50 for Fig. 8 is arbitrary and done only with illustrative purposes. However, 
the behaviour at different multipoles is understandable noting the relative differences among 
different components in the cross-correlations of Fig. 5. This is particularly true for the 
case of the cross-correlation with 2MRS where the curves for the different components are 
remarkably similar at all scales. 

Our main finding is that, although the sterile neutrino auto-correlation component is 
completely dominated by clusters, AGNs and galaxies (see Fig. 3), by taking the cross corre¬ 
lation, we are able to highlight it over the AGNs and galaxies. The contribution from clusters 
of galaxies is very strong also in the cross-correlation case. However, by completely exclud¬ 
ing clusters of galaxies from such an analysis, which is a realistic assumption as eROSITA 
is expected to resolve them all, and reducing the sterile neutrino component analysis only 
to structures up to X h = 10 13 M 0 , we do not loose a significant fraction of sterile 

neutrino decays. In the next section, we will show that a detailed analysis will be able to 
identify the sterile neutrino dark matter over the backgrounds in both cases. 

5 Detectability and projected limits 

In this section, we finally study the sensitivity for sterile neutrino signals using the cross¬ 
correlation angular power spectrum. To this end, we perform a y 2 fit to mock data in different 
energy ranges, with free normalisation of the individual background components and of the 
signal. We take properly into account the covariance of the different components, which is 
usually neglected. 

As already mentioned, we take the soon-to-be-launched eROSITA as reference X-ray 
satellite. Therefore, we adopt an angular resolution = 15" and an (FWHM) energy 
resolution of 138 eV. The angular and energy resolution have a mild dependence on energy 
in the range of interest (see [45] and Footnote 2), which we neglect for this work. The effective 
area is of the order of 500 cm 2 and the field of view (FoV) around 55' diameter, with a strong 
dependence on energy, and we adopt here the energy-dependent values for the grasp as given 
in [45]. The total observation time is assumed to be four years. 
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Figure 7. Contributions to the auto-correlation angular power spectrum for our reference sterile 
neutrino scenario split into different halo mass ranges. Note that for the mass range 10 6 — 10 13 M 0 , 
only the total (1-halo + 2-halo) is shown. Note that our P\i n is obtained for a CDM scenario, therefore 
the contribution below 10 9 M 0 is likely overestimated [67]. However, as clear from the figure, this 
regime has no impact on our conclusions regarding the auto- and cross-correlation power spectra. 




Figure 8. Angular cross-power spectrum for 2MRS (left) and 2MXSC (right) at £ = 50 as a function 
of energy. 

5.1 Statistical method 

The relevant uncertainties that we have to consider in our statistical analysis come from 
the shot noise of the measured X-ray photons, the finite number of measured galaxies, and 
the cosmic variance. The X-ray shot noise affects different X-ray energy ranges differently, 
whereas galaxy shot noise and cosmic variance affect all energy ranges in the same way. 

We perform a sideband analysis fitting the data in three energy ranges simultaneously. 
In this case the cross-correlation angular power spectrum Cj' a becomes a vector with three 
entries, and its variance is described by a 3 x 3 covariance matrix. This matrix is given by 
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(see Appendix A for details) 
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where = 47 r(/) 2 /iV 7 and = 4:ir/Ng are respectively the photon and galaxy count 
shot noise terms, and (I) is the average intensity of the total flux, as shown in Fig. 1, in 
a given energy range. The indices i and j refer to different energy ranges, and 5ij is the 
Kronecker delta. Furthermore, / s ky denotes the fraction of the sky that was observed, and 
Wa = exp(— ct 2 7 2 /2) is the detector window function of a Gaussian point-spread function. 
Therefore, note that the contribution from the photon shot-noise is diagonal, since the noise is 
independent for the various energy ranges, and different in each because the mean intensity is 
different. On the other hand, the galaxy shot noise affects the cross-correlation with photons 
from different energy ranges in the same way. 

We perform our analysis on three reference energies for the sterile neutrino mass: m. s = 
2.0, 7.2 and 18.0 keV. For each of the sterile neutrino masses, we define three energy bands 
that are centred on the line, and are at slightly lower and higher energies. These energy ranges 
are in keV units (low, central, high) = (0.5-0.8, 0.9-1.1, 1.2-1.5), (3.0-3.3, 3.4-3.6, 3.7-4.0) 
and (8.5-8.8, 8.9-9.1, 9.2-9.5). We show the auto- and 2MRS cross-correlation spectra for 
the 0.9-1.1 and 8.9-9.1 keV energy bands in Appendix B. The central band covers most of 
the signal as shown in the left panel of Fig. 8 for 2MRS. In the following, we will discuss 
results only for the 2MRS catalogue for simplicity as the results that one would obtain for 
the 2MXSC catalogue would be less constraining due to the different redshift distribution of 
the latter as explained and showed in the previous sections. 

We treat the instrumental backgrounds in a simple way, by increasing the photon shot- 
noise term, Cjr l , appropriately. We adopt here the factors 5.0, 3.0 and 100 for 2.0, 7.2 and 
18.0 keV sterile neutrinos, respectively. At these energies, the background is dominated by 
Galactic thermal emission, the CXB, and the particle background, respectively [45]. 

The \' 2 function that we use for the sensitivity estimates is given by 


«T - Cg'W) [(Acy^) 2 ]- 1 (CJf - c?fm ), (5.2) 

e=e 0 i,j =i 

where the first term in the first factor denotes the measured cross-correlation spectrum, and 
the second term is the model. The sum is taken over the angular range £ = \£q,£\\ and over 
all three energy bands. Our default angular range, to which we will refer as baseline case, is 
given by £o = 10 2 and l\ = 10 4 . 

The model is a linear combinations of the contributions from sterile neutrinos, unre¬ 
solved AGNs and galaxies, and (resolved and unresolved) cluster emission. Here, we are 
interested in deriving projected upper limits on a sterile neutrino component. We hence 
generate mock data (Asimov data) with the sterile neutrino flux set to zero, and adopt the 
standard Ay 2 method to derive 95% CL upper limits. To this end, we increase the signal 
flux while refitting the other parameters until the y 2 function changes by Ay 2 = 2.71. 


5.2 Projected limits 

The resulting projected limits are shown in Fig. 9. Most interestingly, we find that for our 
baseline analysis eROSITA is sensitive to our reference sterile neutrino scenario, which is 
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projected 95% CL upper limits 



Figure 9. Projected 95% CL upper limits on the sterile neutrino mixing angle as function of the 
sterile neutrino mass. For the baseline results (red solid line) we assume 4 yr of data taking with 
eROSITA, and a sideband analysis of the cross-correlation (CC) angular power spectrum with the 
2MRS catalogue as discussed in the text. This is just on top of the tentative 3.5-keV line reported 
by different groups. For comparison, we also show the limits obtained for 10 yr of data taking (green 
dotted line), and when cross-correlating 4 yr data with a hypothetical perfect model of the dark matter 
decay signal (dashed green line; no shot noise, same window function as for 2MRS and negligible bias). 
We derived projected limits for three reference energies and interpolate otherwise. For comparison, 
we also show regions that are excluded by over- or underproducing the observed DM relic density 
(grey areas), and regions excluded by previous X-ray observations (light brown regions), following 
Refs. [12, 13, 16]. 


motivated by the findings of Refs. [15, 16]. If the claimed 3.5 keV feature is indeed due to 
dark matter decay, eROSITA should be able to see it. 

We discuss now briefly what are the limiting factors for the sensitivity of our baseline 
analysis, and how this could be further improved. In Fig. 9, we also show the limits that 
would be obtained if either the galaxy shot noise or the photon shot noise terms are set to 
zero. These correspond, respectively, to the cases of either perfect knowledge about the dark 
matter distribution (no shot noise, same window function as for 2MRS, and negligible bias; 
note that the latter is close to one for galaxies, see section 3), or 10 yr observation time. We 
find that the galaxy noise term is the most limiting factor for current searches, suggesting 
that improved measurements of the dark matter distribution in the local Universe would be 
of great help when performing cross-correlation searches with eROSITA data. 

Finally, for comparison, we show in Fig. 10 the limits that can be obtained when masking 
out all massive clusters with masses above M 200 x h = 10 13 Mg, for the cross-correlation with 
the hypothetical perfect DM tracer. The limits obtained in this way are only a factor < 2 less 
constraining than what we found with our baseline analysis (with the hypothetical perfect 
DM tracer). This suggest that a good sensitivity to dark matter signals can be obtained even 
if galaxy clusters, with potentially problematic background lines from atomic transitions, 
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Figure 10. Same as Fig. 9, for 4 yr of eROSITA data and cross-correlation (CC) with a perfectly 
modelled dark matter signal. We show the impact of varying the considered angular range (green 
lines), finding that the impact of angular scales £ < 100 is rather minor. We also show how the limits 
change when masking out all halos with masses above M 200 x h = 10 13 M 0 (magenta dash-dotted 
line). The projected limits weaken then by a factor of two, but are then unaffected by the strong 
astrophysical line emission that is associated with the thermal emission from galaxy clusters. 


are excluded from such an analysis as we also showed in Figs. 3, 5 and 7. Note, however, 
that in the current situation where our knowledge of the local dark matter distribution is 
limited, the exclusion of the most massive halos from this analysis is not advisable as can be 
argued looking at our baseline model in Fig. 9: a factor of 2 worsening could compromise the 
ability of testing the 3.5-keV line interpretation. We finally show that including larger angles 
(smaller £) in the analysis would not improve the sensitivity further, whereas the ranges 
l = 10 2 - 10 3 and i = 10 3 - 10 4 are approximately equally important. 

In general, we find that the cross-correlation analysis can significantly improve the 
present constraints on sterile neutrinos below about 10 keV, while above these energies the 
constraints that we obtain are comparable or worse than what can be obtained by pointed 
X-ray observations (e.g., from Refs. [12, 13, 16] as shown in Figs. 9 and 10) as our approach 
becomes dominated by the instrumental background (see also Fig. 12 in Appendix B). 

6 Conclusions 

Motivated by claims of an unidentified 3.56-keV X-ray line found in several galaxies and 
clusters of galaxies, and its interpretation in terms of the decay of sterile neutrino dark mat¬ 
ter [15, 16], we investigated for the first time the possibility of searching for such signatures 
using upcoming full-sky observations of the cosmic X-ray background. We computed both the 
auto- and cross-correlation angular power spectra, where for the latter we adopted a nearby 
all-sky galaxy catalogue, for the sterile neutrino component and astrophysical backgrounds. 
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As astrophysical components that give contributions to the X-ray background, we con¬ 
sidered AGNs, galaxies (powered by X-ray binaries), and thermal bremsstrahlung emission 
from clusters of galaxies. The combinations of these three sources, in particular AGNs, give 
the dominant contribution to the mean background intensity, and the component of the ster¬ 
ile neutrino decays, with inferred parameters from [15, 16], gives only ~1% contribution, even 
at the peak energy around 3.5 keV. 

We found that the auto-correlation angular power spectrum is dominated by the galaxy 
cluster component at almost all angular scales, as they are a rare (more rare than AGNs 
and galaxies) and individually bright sources. AGNs dominate over galaxies and the sterile 
neutrino decays, and the latter gives only a small contribution to the total angular power, 
comparable to that of galaxies. Since eROSITA will resolve most galaxy clusters through their 
thermal emission, the cluster component could be in principle excluded from the analysis. 
Even in that case, however, a large contribution from the unresolved AGNs would hide the 
underlying dark matter component almost completely. 

Since the window function of the spectroscopic 2MASS galaxy catalogue has consider¬ 
able overlap with that of the sterile neutrino decays, taking the cross-correlation with the 
2MASS galaxies efficiently highlight the dark matter component. In fact, we found that the 
cross-correlation power spectrum of sterile neutrino decays is larger than that of both AGNs 
and galaxies. The cluster component is still larger than the sterile neutrino one. However, 
when masking all the clusters that will be resolved with eROSITA, we found that the dark 
matter component (computed from all the halos up to M 200 x h = 10 13 M 0 ) becomes the 
largest in the cross-correlation power spectrum at almost all angular scales. 

Using the above results, we performed a \' 2 analysis in order to estimate the expected 
sensitivity of eROSITA for sterile neutrino searches. We perform a sideband analysis using 
three energy ranges — one centred on the potential line and other two, higher and lower, 
bands to calibrate continuous astrophysical components. We find that when cross-correlating 
eROSITA data with the 2MASS galaxies, the projected 95% CL upper limits after four years 
of observation are cutting into the parameter region of the recently claimed 3.56 keV ster¬ 
ile neutrino decay line [15, 16]. Moreover, we show that the cross-correlation analysis can 
significantly improve the present constraints on sterile neutrinos below about 10 keV, while 
above these energies our approach becomes dominated by the instrumental background and 
the corresponding constraints are not as competitive with respect to pointed X-ray obser¬ 
vations. We also find that the masking of galaxy clusters does not improve the constraints, 
unless a better knowledge of the local dark matter distribution is achieved. In fact, the main 
limitation comes from the shot noise associated with the finite number of galaxies in the 
2MASS catalogue. Indeed, we showed that when using shot-noise free tracers for dark mat¬ 
ter, i.e., having a perfect knowledge of the local dark matter distribution with no bias, the 
projected limits improve significantly, allowing eROSITA to put stringent constraints on the 
possible sterile neutrino signal. Hence, in order to ensure that the potential of eROSITA and 
other full-sky X-ray surveys is fully realised for dark matter searches, it is critical to obtain 
a detailed description of the dark matter distribution at low redshift, at galaxy and galaxy 
cluster scales, without the shot-noise limitations of the 2MASS catalogue. To what extent 
this can be achieved by a combination of optical and X-ray observations (e.g., Ref. [32]), and 
of gravitational lensing [29, 33] will be discussed in a forthcoming publication [74]. 
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Note added. In the final stages of our work we became aware of another group exploring 
a similar approach [75]. 
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A Treatment of correlated statistical uncertainties 


We discuss briefly how the covariance matrix used in our y 2 analysis follows from first princi¬ 
ples. For a discussion about the angular power spectrum in light of finite counts see also [76]. 

Our analysis is based on the cross correlation angular power spectrum that is obtained 
from the measurements in some X-ray energy band i and the galaxy catalogue g (mostly 
2MRS), and given by 


— ( a Jmi( a em) ) — 
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Here, aj mi («| m ) are the modes of a decomposition of the measured X-ray sky map (galaxy 
catalogue) into spherical harmonics, and (...) denotes the usual average over m. 

We are interested in the variations of Cjf over many measurements. In general, it holds 
that ( a lmi)r — ( a lm) r = 0, where (. ,.) r is the average over cosmological realisations, and 
photon and galaxy shot noise. Furthermore, we find that 
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where IV 7i i is the number of photons in energy bin i, and (/)$ the corresponding full-sky 
averaged flux. The noise comes either from the finite number of photons in X-ray observations 
or the finite number of galaxies in the catalogue, and the cosmic variance contribution is a 
consequence of the stochasticity of structure formation. We assumed that the morphology of 
the signal changes only weakly between the relevant energy bins. 

The covariance of the cross correlation angular power spectrum is then obtained as 
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Figure 11. Left. Auto-correlation angular power spectrum of CXB due to sterile neutrino decays for 
m Va = 2 keV with a mixing angle corresponding to the constraint given by our baseline analysis as in 
Fig. 9, (unresolved) AGNs, (unresolved) galaxies, and (resolved and unresolved) clusters of galaxies 
in the 0.9 — 1.1 keV energy band. For AGNs and galaxies, both the correlation term (solid) and 
total including Poisson term (dotted) are shown. Right. The same but for m Vs = 18 keV in the 
8.9 — 9.1 keV energy band. We show both plots with the same scale for comparison, but note that 
the scale is different with respect to the left panel of Fig. 3. 


The adopted Eq. (5.1) is then obtained after a rescaling of the photon power spectrum 
that inverts the (here energy-independent) effect of the finite angular resolution, e.g., —> 

cl/wj. 

B Auto- and cross-correlation power spectra 

We show in Figs. 11 and 12 the auto-correlation and 2MRS cross-correlation power spectra in 
the 0.9 — 1.1 and 8.9 — 9.1 keV energy ranges corresponding to the sterile neutrino masses of 
m Us = 2 and 18 keV, respectively. In both cases, the sterile neutrino mixing angle is chosen 
to be just as from the constraints of our baseline analysis shown in Fig. 9. Our constraining 
power decrease dramatically moving toward higher energies, as clear both from Fig. 9 and 
from the error bars in the left panel of Fig. 12, due to the larger instrumental background. 
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Figure 12. Cross correlation with 2MRS. We overlay to the cluster component the error bars, as 
A Ci/yfl/2, estimated from the diagonal part of the covariance matrix of Eq. 5.1. Left. Cross¬ 
correlation angular power spectrum of CXB due to sterile neutrino decays for m Us = 2 keV with a 
mixing angle corresponding to the constraint given by our baseline analysis as in Fig. 9, (unresolved) 
AGNs, (unresolved) galaxies, and (resolved and unresolved) clusters of galaxies in the 0.9 — 1.1 keV 
energy band. For sterile neutrinos, we also show the case when the integration upper mass limit is 
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